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1. Introduction 

Starting from the one-loop level, calculations of the squared amplitudes in Quantum Field 
Theory bring us to consider rc-point tensor integrals of rank R, (n,R) -integrals 

J in^ n ;' =1 c7 ' ( } 

where d = 4 — 2e and denominators cj have indices Vj and chords qj 

Cj = (k-qj) 2 -mj + is. (1.2) 

There are many strategies how to solve such tensor integrals, see e.g. [1-4] and references therein.. 

Our procedure and notations are described in [5, 6] and in references quoted therein. There 
are two main steps which are applied to (1.1). First, higher dimensional scalar integrals are used. 
This is Davydychev's decomposition [7], e.g. for a R = 2 tensor with a general number of external 
legs n 

# V = f d kn v f[c- X = £ rfffjvy^-l^^ 1 - (1-3) 

J r=\ ij=l Z 

We use the notation 

[d+]'=d + 2l. (1.4) 

Second, recursive relations as derived in [8, 9] bring the expression into a simpler form, e.g. 
explicitly for pentagons 



U5 s=l,^i \>5 \h 



As may be seen, there is a decrease in both dimension and indices. Finally, recursively, we go 
down to basic scalar integrals in higher dimension - but at the price of appearance of inverse Gram 
determinants () 5 in (1.5). It is known that inverse Gram determinants can cause trouble in realistic 
physical applications, because there are kinematical regions where they can be small or even vanish. 

Here, the main aspect of the new approach comes into action: After insertion of (1.5) into 
(1.3), the chords qi,qj are contracted with external momenta (expressed as simple differences of 
chords q a ,qb)- One observes the appearance of "auxiliary vectors" 

a v = E*}|r. * = o,..,5. (1.6) 

j=\ Us 

Contractions of the vectors with a chord q a immediately eliminate the unwanted inverse of () 5 . At 
first glance this is really surprising and we want to stress this property in particular. If we would 
arrange cancellations of inverse Gram determinants () 5 in (1.3) in a standard way, much more effort 
is needed [10, 11, 5]. 
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To be more precise, let us give a manifest example, taken from subsection 2.1: 

"- 1 ('•) 1 

(?«•&) = 7T = o^- 5 ^). * = !,.■■, 5- (1.7) 

7=1 U » / 

The 2 vectors, when contracted with a chord in (1.7), give just Kronecker 5-functions. The inverse 
Gram determinant () 5 , present in intermediate steps, disappears automatically. 

Let us account now the main advantages of the new approach where chords are contracted 
with <2 vectors (signed minors): 

1. Cancellation of inverse Gram determinants () 5 . 

Explained as above - occasionally also cancellation by subtraction. 

2. Elimination of the metric tensor g^ v . 

The integral (1.3) in its original form includes the g^ v metric tensor. 
The identity (valid in d = 4 dimensions) 

1 5 (•) 

r=If^ ; V (1-8) 

introduces the inverse of () 5 again. So, as argued before, it should cancel out. In fact, e.g. 
for pentagons (n = 5), 

ir=i<iU)v ij if;/-\ 8 ^i [ t\ (i.9) 

and the last term in (1.5) cancels against the g^ v term in (1.9). After the g^ v cancellation and 
the contractions with chords, the remaining inverse of Q 5 , which was produced in interme- 
diate steps, disappears as well. 

3. Presence of Kronecker5-f unctions and simple kinematical factors Y . 

This is probably the main feature of the approach. Due to relations like (1.7) and 

n-\ (?) ! 

-Go) = Y,^^i)ir = -0^- Y nn), (1-10) 
7=1 U « Z 

where Yy = —(qi — qj) 2 -\-mf + mj are simple kinematical functions, the expressions for 
contracted tensor integrals presented in the next section take a simple form. Take as an 
example the (2.33) from the next section. For given indices a,b and t, most of the terms 
vanish. 

The features are independent of multiplicity (i.e. the number of external legs) and rank of the tensor 
integrals under consideration. 

2. The contractions 

We consider now 5-point functions and eliminate the inverse 5-point Gram determinants () 5 . 
Inverse 4-point Gram determinants Q g will be treated separately. The explicit formulae cover 
tensor integrals until rank R = 3. 
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2.1 The vector integral: /j 

The vector integral may be written in terms of the scalar 5-point function E and 4-point scalar 
functions /|, where s indicates the scratched line of the 5-point topology: 

$=E$-£W, E = ^-£( S )P 4 . (2.1) 

For q„ = 0, a = 1, . . . ,n — 1, s = 1, . . .n, the multiplication with chords introduces (1.7) and 
(1.10), and the final expression for a contraction of the vector integral /r with an arbitrary chord 
q a ^ becomes 

5 

(q a -h) =E(q a -Q )-Y J ll{qa-Q s ) 

s=\ 

= -\e (Y a5 - Y 55 ) - - £ II (8 as - 8 5s ) . (2.2) 

Here we have applied (1.7) and (1.10) for n = 5. In fact no inverse Gram determinants have 
appeared. 

2.2 The tensor integrals of rank R = 2: 

From equation (3.5) of [5] one gets 

with 

There are two kinds of contractions to be considered: with the metric tensor or with two chords. In 
principle, one may also meet in chiral theories contractions of the form 

Q 1 Zapvv h ■ ( / -- ) ) 

Proper organization of the matrix element evaluation will exclude the appearance of such terms 
since cross sections are either helicity and spin summed and thus pure scalar quantities or written 
in terms of helicity states. In the latter case, the scalar factors with the loop functions are also free 
of 75, and the e tensor might influence the matrix element basis. Nevertheless, we add one rank 
R = 2 example in order to demonstrate that this kind of sums may also be treated. 

2.2.1 Contractions with chords: q a \iqb\^ y 
By contraction, (2.3) rewrites into 

5 

q a ^qbyis V = {q a -h){qb • Go) - £ " 7 4) {lb • Q*)- ( 2 -6) 
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Here no inverse () 5 occurs anymore. Inserting the sums and E a ' rt from [6], {q a -/J) can be 
written as 



where we introduced the abbreviation 



(Y a5 - YssYa + - 8 5t )lf + (S as - 8 5s )R s \ ,(2.7) 



1 

5 



'4-1 



f=l 



r.vr 
'3 



(° S ) 
\0s)i 



t=\ 



oj; 4 



t=i 



rSt 
y 3 



(2.8) 



The 2 nd representation of R s in (2.8) can be used if the four-point Gram determinant Q 5 is small 
or even vanishes. 

2.2.2 Contractions with the metric tensor: g^ v /^ v 

The second scalar which can be constructed from the tensor of rank 2 is g^ v I^ V ■ Due to (2.3) 
we have to evaluate the following scalar products: 



(Go -Go) 

(Go-G,) 

(GS-GO 
(Gt • GO 





(o) 



5 -+Y 55 



() 5 

S S 5, 



5 - - au 



In this case the terms with J- cancel and, not surprisingly, the result finally is 



g^5 



2.2.3 Contractions with the antisymmetric tensor: q a,1 q bv s^ va pl^ 



(2.9) 



(2.10) 



A reduction of pseudoscalar contractions leads to expressions of the type 

P[C)=qPqie tlvpK £ qfq]Qj. 



(2.11) 



This may be evaluated as follows: 



P[C] = qPq X b W g\, g\, £ qf'qfa 



01 



t,WV/S »n,n=l\ n /5 ij=l 



4 4 

a+i+fe+w=10 



(fall £flV P Xlk1m 



!,n=l 



£ {<irn)(<in-qj)Ci 4 



4 





2 det [91, ? 2 , 93, 94] £ S(a,b,k,m) £ f m> ) I/„[C], 

5 fc,m=l /,«=! V/5 W 5 



(2.12) 



a+fo+£+m=10 
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where S(a,b,k,m) gives the sign of permutations needed to place the indices in increasing order. 
Further we introduced the abbreviation 



E /n[c] = £ {qvqi){qn-qj)Cij, 

the calculation of which has to be performed for specific cases. 
Let us choose for demonstrational purposes Cn = ) . 



(2.13) 



Zl n {C] = -{qvq n ) 



-() 5 (Si s -8 5s )(8 ns -8 5s ) 



(2.14) 



i(qrqn)O 5 -l0 5 s is s ™ s ^5, 
i(qrqn)O 5 -l0 5 s = 5 - 



Next we have to evaluate the sum (;) 5(^)5 ' Taking the first term in (2.14), to begin 

with, we have 



* 1=1 
Qsfk 



1 I 5 



4 V m 



and for s 7^ 5 we thus have 



'n=l 



1 1 r\n 



in 



■E/n[C] 



05 

4 



4 Vww 



Similarly we have for s = 5 



t.c 



•£„[c] = % 



111 



5 /«=1 



Due to 



1=1 



we get for 5 = 5 the same result as in (2.16). 

Thus the total result reads for the specific choice of Cy): 



(2.15) 



(2.16) 



(2.17) 



(2.18) 



is 



. =det[gi, 92,^3,94] XI S(a,b,k, 

3 s J 5 jfc,m=l 

a+i>+fc+m=10 



/v.V 
H/.V 



(2.19) 
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In the summation of (2. 19) all indices a,b,k, m must be different, i.e. the sum contains only 2 terms, 
with k and m exchanged. Since these terms come with opposite sign and (2.16) is symmetric in k 
and m, the final result for the present example vanishes, i.e. 

4 



0. 



2.3 The tensors ofrank/? = 3: /; 

From eq.(3.5) of [5] one gets 



(2.20) 



s=l 



where according to (3.12) in [5] it is 



k!r ^ ti.v ti.v y 4 ' 



(=1 



y 3 



«=i 



(2.21) 



(2.22) 



2.3.1 Contractions with chords: qa^qbvq c x^ Vl 

The contraction of (2.21) with ^ is trivial due to (1.7) and (1. 10). Further, the term q a iiqbv^ 
is known from subsection 2.2.1. The only new object to be investigated is therefore 

qarfbvir* = (?« ■ I'D (?* ■ 2 s o) " £ (q a ■ /?) (q b ■ G?) - jk- (qa ■ &)(?* ■ Qs)^ ■ (2.23) 

W5 



r-1 



New sums may be expressed by those derived in [6] 

(qt-Qo) 



1 y 2,s 
J 5 



and 



fa I?) 



1 



©5 



v3,.vf jrf 
y 3 



(q b -Q s t ] 



Ly\,Stu jStu 
Z-'a n 

»=1 



1 

0: 



yl,St 



H=l 



where we introduced as a further abbreviation 

1 



©5 
1 

\0stJi 



jStU 

l 2 



U '/ 5 «=1 



0s? 
Owf 



JStU 



(2.24) 
(2.25) 



- Fss)/? + £ («« - Ss u )I^ + (8 as - 8 5s )R st + (S at - 8 5t )R ts \ ,(2.26) 



(2.27) 
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From 



Am 



(8 bt - 8 5t ) 



(8 bs - 8 5s ) 



(2.28) 
(2.29) 



multiplying (2.7) and (2.26), respectively, we see that the (8 bs — 55 S )-terms in (2.29) combine with 
the first terms in (2.7) and (2.26) to yield the contribution 



1 



(Y a5 -Y 55 )(8 bs -8 5s )R s . 



(2.30) 



Collecting all contributions, the following relation is useful in order to to eliminate some t4 — terms: 



0, 



( s ) 



( S ) 5 ( S ) 

7^r 7 3 - L T^r 72 • 



(2.31) 



KsJS \s)5 u=\ \s)5 

The result for (2.23) is explicitly symmetric in the indices a,b: 

1 



jHV,s i,~ 2 \/s s \ O5 i[d+]- s 

q uv ,qbvl\ = q ai iqbvh - 7 \8 as - 8 5s ) (o bs - b 5s ) 777-/4 , 



(2.32) 



and the first term on the right hand side is 



qa\iqbv L A 



pv,s 



^(Y a5 -Y 55 )(Y b5 -Y 55 )-I s 4 
+ 1 {8 as - 8 5s )(8 hs - 8 5s ) 
+ 7 [(Sas ~ 8 5s )(Y b5 -Y 55 ) + (F a5 - Y 55 )(8 hs - 8 5s )]R s 



1 



£ [(Sas - 8 5s )(8 bt - 8 5t ) + (8 bs - 8 5s )(8 at - 8 5t )]R s 



£ [(Sat ~ 8 5t )(Y b5 - Y 55 ) + (Y a5 - Y 55 )(8 bt - 8 5t )]L 



1 



r=l 
5 



1 



Z(S at -S 5t )(S bt -S 5t )R ts + - £(S bt -8 5t )(8 au -8 5u )I?, tu . (2.33) 



^ t=l f,U=l 

It is worth mentioning that the second term in (2.32) may be written as 

5 



(U 



R s . 



(2.34) 



The R s , R st , £j =1 (j) 5 ^ f and the scalar integrals in (2.33) are independent of the indices a,b,c 
and can be considered as buildung blocks. The summation over s,t,u is simplified by the explicit 
appearance of the Kronecker 5-functions, i.e. all contractions are just given in terms of different 
combinations of the buildung blocks. 
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There still occur contributions with factor vi- in (2.32) and in (2.33). At first we observe that 

V.J 5 

this contribution is indeed quite "localized": in general one will be able to avoid the small Q to 
appear on line s = 5 by choosing the integration momenta properly. Let say Q°) is small. Then 
the critical contribution will only occur if a = b = c = sq, i.e. only in a single contraction. 

Nevertheless, if one wants to avoid the -Jr- completely, one can express the corresponding 

W5 

contribution for the sum of these terms in (2.32) and (2.33) in terms of higher dimensional integrals, 
making again use of (2.8). Let us call the sum 7|: 



Ji = 




■o 5 i [ r ] ' s + 



5 t=\ 



R s 



5 (d-2)(d-i)ir r - s - Qit u +L (^-2)4^+1 



(2.35) 



The higher dimensional integrals can then be evaluated as described in [5]. 

Let us call the second term on the right hand side of (2.21) C$. abc . It becomes, when contracted 
with momenta q a , q b and q c : 



Cs.abc = - qafiqbv<lcX 2-4 &s 
s=l 

5 1 



£ ^a\iqbv^ ,s ■ ~ (Scs - S 5s 



(2.36) 



The last term again contains (2.32). Introducing Y a = Y a s — F55 and making use of (2.35), we may 
collect terms and get 

1 



'5,abc 



{j 5 4 - B ab 8 ac Jl + Y fl Y, (7 4 5 - II) + Y a (/f 5 + if - if) + Y, (if + if - If) 



1 rab5 I racj , j, 
-\-i 2 -r h 



ac5 



T bc5 



jcibc 



Y a (R 5 + 8 bc R b )-Y b (R 5 + 8 ac R a ) 



. R 5a_ R 5b_ R 5c + g 



'ab 



R 



a5 



R aC )+8 a , 



R a5 R ab 



+ 8 bc {R b5 -R h " 



(2.37) 



2.3.2 Contractions with the metric tensor: gp.vq.axl5 



fivX 



mvX 



In order to calculate g^vl^ v '\ we need g^vl^ and thus the further scalar products, see (2.22): 



(G5 • OS) 
(Q s -Qs) 
(Ql ■ Q s o) 



m ■ at) 



1 



2C) 5 

105 
2 5 ' 
1 

W 5 



W5 



(G? ■ G?) = 0, 



+2 



->v5 



1 

+ 2 755 ' 



8 t5 + 



8 t5 + 



%5 



Js5 



(2.38) 
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which yields 
and finally 



Sjiv/T' • If I 8 s5 R 5 (2.39) 



*Mv?«a4"* = "X [(Ya5-Y 55 )E+I a 4 -II] - \ [(Y a5 -Y 55 )I 5 4 +If -R 5 ] . (2.40) 

It is remarkable that (2.39) is trivial again for s ^ 5. For s = 5, however, the standard cancellation 
of propagators does not work and for this case (2.39) is indeed a useful result. 

3. The OLEC package - first numerical results 

The relations introduced in the previous sections have been implemented in Mathematica. A 
corresponding file OLECv0.9.m can be found at [12]. Certainly, computer algebra programs like 
Mathematica are not the optimal and efficient tools for ultimate numerical applications, they are 
slow. However, there are some reasons to start with Mathematica. First, standard tensor reductions 
for five and six point functions discussed in [11] have been also implemented in the Mathematica 
package hexagon.m [13, 14]. The new method of direct contractions of tensors as introduced 
in [6, 5] and presented here can be compared with hexagon.m. Second, if basic scalar integrals 
are known analytically, then the complete reductions are automatically given also in an analytic 
form. This approach could be used for instance to investigate infrared structure of amplitudes 
[15-17, 11]. Of course, for numerical applications one will prefer packages written in Fortran or 
C++. For efficient reductions of five-point tensor integrals, see V. Yundin's C++ package PJFry 
[18-20]). 

For the contracted 5-point functions with which we are dealing here, the following notation 
is used: CEx (chords) , where "C" stands for "contracted" and "E" stands traditionally for pen- 
tagons. The argument "chords" represents a list of indices, and their length depends on the rank "x" 
of the considered object. In this way, CE\{a) is equivalent to (2.2), CE2(a,b) is equivalent to (2.6) 
and CE3(a,b,c) is equivalent to appropriate relations given in subsection 2.2.1. For the objects 
contracted with the metric tensor g^y we have two objects up to rank 3: CE2g and CE3g(a). They 
correspond to (2.10) and (2.40), respectively. 

Now, we can compare these objects with corresponding objects which come from standard 
reductions. Below we will use the notation of LoopTools/FF [21, 22] for coefficients of corre- 
sponding tensorial objects EOi: 

CE\(a) = qfr j J^$j 11 E0i[eei,iwj] j , (3.1) 

CE2(a,b) = q%'q£ 2 I £ q^q^EOileeij ,invs\ + g^^ 2 E0i[ee00,invs] > , (3.2) 
CE3(a,b,c) = cfi'qfq? 3 I £ q'^q^q^EOileeijkJnvs] 

+ E (Sm^m +g^ 3 q'n 2 +gn 2fi3 q' fll )E0i[ee00i,invs] 1 . (3.3) 
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For contractions with the metric tensor: 

CElg = E0i[ee00l, invs] , (3.4) 

CE3g(a) = g^q% J £ q l lii q } fl2 q J ^EOi[eeijk,invs] 
{{iJM=l 

+ L (Smm^ + Suited + Stm^in ) E0i[ee00i, invs] j . (3.5) 

Numerical results are shown in the file Contract s_examples . nb. They are compared for 
various chord indices using the Looptools/FF and OneLoop [23] packages. 

Below we tabularize an example of results extracted from Contracts_examples . nb for 
kinematical points defined by: 
p\ = p 2 2 = 0; p\ = p\ = 49/256; pj = 9/100; 
*i2 =4; ^23 = -1/5; S34 = 1/5; s 45 = 3/10; 515 = -1/2; 
m\ = m\ = m\ = 49/256; m\ = m\ = 81/1600. 

In: 

Do[Print[i, ") LoopTools: ", RanklLoopTools [ i ] , ", Hexagon: ", 

Ranklhexagon [i] , ", Contracts: ", CE1 [i] ] , {i, 1, 5}] 
Out : 

1) LoopTools: -19.1597+12.6772 I, Hexagon: -19.1597+12.6772 I, 
Contracts: -19.1597+12.6772 I 

2) LoopTools: -80.4038+91.9442 I, Hexagon: -80.4038+91.9442 I, 
Contracts: -80.4038+91.9442 I 

3) LoopTools: -19.1525+12.4364 I, Hexagon: -19.1525+12.4364 I, 
Contracts: -19.1525+12.4364 I 

4) LoopTools: -17.496+10.67060 I, Hexagon: -17.496+10.67060 I, 
Contracts: -17.496+10.67060 I 

5) LoopTools: 0, Hexagon: 0. +0. I, Contracts: 



In addition, comparisons concerning the speed of calculation have been made. As expected, the 
Mathematica package OLEC.m v.0.9 is at least one order of magnitude faster than hexagon.m but 
still slower than the Fortran package LoopTools/FF. 

Thus, a second package OLEC v.0.9 written in C++ has been prepared for numerical studies. 
It should be assumed as the skeleton for the first version of the ultimate package. It has to be further 
optimized and additional refinements for specific kinematical points are needed, similar to those 
which were implemented for the standard reductions within the PJFRY package. 

The OLEC package can be downloaded from [12]. It is linked with both the LoopTools/FF 
and the OneLoop libraries. Our first package, where the contractions replace the explicit tensor 
reductions is at the moment comparable in speed to numerical results of 5-point tensor reductions 
implemented in LoopTools/FF. This has been tested for specific rank 2 tensors in (3.2) and tensors 
of rank 3 in (3.1). Details can be found in the module example_2 of the Makefile in [12] where 
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CE2(l,2) and C2s3(l,2,3) can be calculated and compared with corresponding LoopTools objects 
using different flags, for users convenience. The result is not too surprizing. Evidently, for higher 
tensor ranks the advantage of contractions compared to rank reductions, introducing multiple sums, 
should become more pronounced. 

4. Conclusions and outlook 

A new method of evaluation of contracted tensor integrals has been presented for five-point 
one loop integrals. So far individual objects like CEl, CE2 and CE3 have been resolved and 
implemented into first codes. What remains is to apply it to the calculation of whole amplitudes 
for some physical processes. 

In the near future the first version 1 .0 of the complete code in C++ and Fortran 90 is planned 
to be released. Here we have shown the backbone of the code with first numerical estimations. 
In addition, higher than rank 3 tensor structures for 5-point functions and tensorial integrals with 
more than five external legs are in preparation. 

Acknowledgments 

We would like to thank Andreas van Hameren, Thomas Hahn and Krzysztof Kajda for use- 
ful discussions and help. Work supported by European Initial Training Network LHCPHENOnet 
PITN-GA-2010-264564 and by Polish Ministry of Science and Higher Education from budget for 
science for 2010-2013 under grant number N N202 102638. 

References 

[1] Z. Bern et at, The NLO multileg working group: Summary report, 2008. Presented at 5th Les 
Houches Workshop on Physics at TeV Colliders, Les Houches, France, 11-29 June 2007. 
[arXiv:0803.0494]. 

[2] J. Andersen et at, The SM and NLO Multileg Working Group: Summary report, 2010. Proceedings of 
the 6th Workshop on Physics at TeV colliders, 8-26 June 2009, Les Houches, France. Dedicated to 
Thomas Binoth. [arXiv:1003.1241]. 

[3] J. Alcaraz Maestre et ai, The SM and NLO Multileg and SM MC Working Groups: Summary Report, 
2012. Report of the SM and NLO Multileg and SM MC Working Groups for the Workshop 'Physics 
at TeV Colliders', Les Houches, France, 31 May - 8 June, 2011. See [arXiv: 1203. 1488] and 
[arXiv:1203.6803]. 

[4] J. Bliimlein, S. Moch, and Tord Riemann (eds.), Loops and Legs in Quantum Field Theory, 

PoS(LL2012) (2012). Proceedings of the 1 1th DESY Workshop on Elementary Particle Theory, April 
15 - 20, 2012, Wernigerode, Germany. [PoS(LL2012)]. 

[5] J. Fleischer and T. Riemann, Complete algebraic reduction of one-loop tensor Feynman integrals, 
Phys. Rev. D83 (201 1) 073004, [arXiv : 1 9 . 4 4 3 6]. 

[6] J. Fleischer and T. Riemann, Calculating contracted tensor Feynman integrals, Phys. Lett. B701 
(2011) 646-653, [arXiv: 1104 .4067]. 



12 



Contraction ofN-point tensor integrals 



Jochem Fleischer 



[7] A. I. Davydychev, A simple formula for reducing Feynman diagrams to scalar integrals, Phys. Lett. 
B263 (1991) 107-111. doi:10.1016/0370-2693(91)91715-8. 

[8] O. Tarasov, Connection between Feynman integrals having different values of the space-time 
dimension, Phys.Rev. D54 (1996) 6479-6490, [hep-th/ 9 60 6018]. 

[9] J. Fleischer, F. Jegerlehner, and O. Tarasov, Algebraic reduction of one-loop Feynman graph 
amplitudes, Nucl. Phys. 8566(2000) 423^40, [hep-ph/9907327]. 

[10] J. Fleischer, Application of Mellin-Barnes representations to the calculation of massive five-point 
functions in bhabha scattering, 2007. Talk held at Conference on Frontiers in Perturbative Quantum 
Field Theory, 14-16 June 2007, ZiF, Bielefeld, Germany. Transparencies: 
fleischer-bielefeld-2007-talk.pdf. 

[11] T. Diakonidis, J. Fleischer, J. Gluza, K. Kajda, T. Riemann, and J. Tausk, A complete reduction of 
one-loop tensor 5- and 6-point integrals, Phys. Rev. D80 (2009) 036003, [arXiv : 8 12 . 2 1 3 4]. 

[12] J. Fleischer, J. Gluza, M. Gluza, and T. Riemann, OLEC project; C++ and Mathematica packages 
OLEC v.0.9 available at Silesian Univ., Katowice, webpage http://prac.us.edu.pl/~gluza/olec/. 

[13] J. Fleischer, J. Gluza, K. Kajda, T. Riemann, and V. Yundin, Hexagon project; Mathematica package 
hexagon. m v. 1.0 available at Silesian Univ., Katowice, from 
http://prac.us.edu.pl/~gluza/hexagon/packages/hexagon.rn. 

[14] J. Gluza, K. Kajda, T. Riemann, and V. Yundin, New results for loop integrals: AMBRE, CSectors, 
hexagon, PoS ACAT08 (2008) 1 24, [0902.4830]. 

[15] J. Gluza and T. Riemann, New results for 5 -point functions, eConf C0705302 (2007) LOOP01, 

[0712 .2969]. 

[16] J. Gluza and T. Riemann, A New treatment of mixed virtual and real IR- singularities, PoS 
RADCOR2007 (2007) 007, [ 8 1 . 4 2 2 8 ] . 

[17] T. Diakonidis, J. Fleischer, J. Gluza, K. Kajda, T. Riemann, and J. Tausk, On the tensor reduction of 
one-loop pentagons and hexagons, Nucl. Phys. Proc. Suppl. 183 (2008) 109-115, 

[arXiv: 0807 .2984]. 

[18] J. Fleischer, T. Riemann, and V. Yundin, PJFry: A C++ package for tensor reduction of one-loop 
Feynman integrals, . In: [3], preprint DESY 11-25 (2011), 
http://www-library.desy.de/cgi-bin/showprep.pl7desyl 1-252. 

[19] J. Fleischer, T. Riemann, and V. Yundin, New developments in PJFry, 1210.4095. 

[20] V. Yundin, Massive loop corrections for collider physics, PhD thesis, Humboldt-Universitat zu Berlin, 
2012, see http://edoc.hu-berlin.de/docviews/abstract.php?id=39163 and 
urn:nbn:de:kobv:l 1-100199626. 

[21] T. Hahn and M. Perez- Victoria, Automatized one loop calculations in four-dimensions and 
D-dimensions, Comput.Phys.Commun. 118(1999) 153-165, [hep-ph/9807565]. 

[22] T. Hahn, Feynman Diagram Calculations with FeynArts, FormCalc, and LoopTools, PoS ACAT2010 
(2010)078, [100 6.2231]. 

[23] A. van Hameren, Oneloop: For the evaluation of one-loop scalar functions, Comput. Phys. Commun. 
182 (201 1) 2427-2438, [arXiv =10 07.471 6]. 



13 



